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Abstract 

Energy levels of neutral atoms have been re-examined by applying an alterna- 
tive perturbative scheme in solving the Schrodinger equation for the Yukawa 
potential model with a modified screening parameter. The predicted shell 
binding energies are found to be quite accurate over the entire range of the 
atomic number Z up to 84 and compare very well with those obtained within 
the framework of hyper-virial-Pade scheme and the method of shifted large- 
N expansion. It is observed that the new perturbative method may also be 
applied to the other areas of atomic physics. 
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I. INTRODUCTION 



In recent years the energy levels of neural atoms have been studied by several analytic 
methods [1-8] in which it is assumed that the screened potential of the atom may be of 
static screened Coulomb (SSC) which is well represented by Yukawa form: 

^(r) = -(^)exp(-5r), (1) 

with A = aZe"^, where a = (137.037)"^ is the fine-structure constant and Z is the atomic 
number. This form is often used for the description of the energy levels of light to heavy 
neutral atoms [7]. It is known that SSC potential yields reasonable results only for the 
innermost states when Z is large. However, for the outermost and middle atomic states, 
it gives rather poor results. Although the bound state energies for the SSC potential with 
Z — 1 have also been studied. The screening parameter 5 is chosen to be 

S = SoZ'f', (2) 

corresponding to the Z-dependence of the reciprocal of the Thomas-Fermi radius of the 
atom. However, these analytic works, perturbative as well as nonperturbative, fail to yield 
accurate shell binding energies for light atoms, particularly in the range Z <9. Subsequently, 
it has been pointed out by Refs.[5,6] that the major source of errors perhaps lies in the wrong 
choice of the Z-dependence of the screening parameter. Invoking Fermi- Amaldi correction 
[9] in the context of Ecker-Weizel approximation (EWA) method [10], Dutt and Varshni 
[5,6,7] have suggested a modified form 

S^SoZ^/^l-l/zf\ (3) 

with 6q = 0.98. Clearly, when Z = 1, 6 vanishes and the potential in (1) becomes the 
Coulomb potential as it should be. Correctness of the choice of the modified screening 
parameter has been further justified by the recent work of Lai and Madan [8]. They have 
shown that the hypervirial-Pade scheme which failed to reproduce correct shell binding 
energies for light atoms using the screening parameter given in (2) [4] , yields very accurate 
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energy eigenvalues using the modified screening coefficient in (3) [8]. However, one of the 
short-comings of the hypervirial-Pade technique is that it involves elaborate computational 
time and effort for each numerical prediction. Lai and Madan [8] have to consider upto eleven 
terms in the perturbation series for the energy eigenvalues in order to ensure the convergence 
of the Pade approximant -^(A^, M). Furthermore, application of this method becomes quite 
restricted due to nonavailability of compact analytic expressions for the bound-state energies, 
eigenfunctions and normalization constants. 

On the other hand, Dutt and Varshni [7] have investigated the bound-states of neutral 
atoms using the large-A^ expansion method which has been claimed to be very powerful 
for solving potential problems in nonrelativistic quantum mechanics. This technique also 
requires an approximate treatment and computational time as well. 

In this paper, we investigate the bound-state properties of SSC potential using a new 
perturbative formalism [11] which has been claimed to be very powerful for solving the 
Schrodinger equation to obtain the bound-state energies as well as the wave functions in 
Yukawa or SSC potential problem [11,12] in both bound and continuum regions. This novel 
treatment is based on the decomposition of the radial Schrodinger equation into two pieces 
having an exactly solvable part with an addiitional piece leading to either a closed analytical 
solution or approximate treatment depending on the nature of the perturbed potential. 

It seems then logical and meaningful to probe whether the range of applicability of this 
novel perturbation treatment may be widened. As a first attempt, we have shown recently 
that the method adequately explains the spectrum of hydrogen-like atoms A = Z = 1 and 
also light and heavy atoms [13]. With a view to make further applications to problems of 
atomic physics, we compute here the shell binding energies of light to heavy neutral atoms. 
The relevant steps of the perturbation scheme are to obtain analytical expressions for the 
bound-state energy levels and corresponding normalized eignf unctions. 

The contents of this paper is as follows. In section II we breifiy outhne the method with all 
necessary formulae to perform the current calculations. In section 111 wc apply the approach 
to the Schrodinger equation with SSC potential and present the results obtained analytically 
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and numerically for the bound-state energy values upto third perturbation energy shift. 
Finally, in section IV we give our concluding remarks. 



II. THE METHOD 



For the consideration of spherically symmetric potentials, the corresponding Schrodinger 
equation, in the bound state domain, for the radial wave function reads 



2m ijjn (r) 



V{r) - K, 



(4) 



with 



Voir) + 



+ 11 



2m 



+ AV{r), 



(5) 



where AV{r) is a perturbing potential and ipni^) — Xnir)un{r) is the full radial wave func- 
tion, in which Xn{r) is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas Un{r) is a moderating wave function corresponding to the perturbing po- 
tential. Following the prescription of Refs. [11,12], we may rewrite (4) in the form: 



V{r) - En. 



(6) 



2m \Xnir) Un{r) Xn{r)un{r) ^ 

The logarithmic derivatives of the unperturbed Xn{r) and perturbed Un{r) wave functions 
are given by 



^Xn{r) 



and AW„ 



^ Kir) 
^2mUn{r)' 



(7) 



which leads to 



W^ir) - 



h 



-W^ir) 



Voir) + 



iii + 1) 

2m 



- £r. 



2mxnir) ' V2m 

where e„ is the eigenvalue for the exactly solvable potential of interest, and 

h 



(8) 



f u'nir) ^ r, Xnir)Kiry 



2m V u„(r) 



Xn[r)Un[r) 



AW^ir) - ^=A<(r) + 2Wnir)AWnir) = AT/(r) - Asn, 



V2 



m 



(9) 



in which A£„ = Ej:^^ + Ej^^ + E^^ -\ is the correction term to the energy due to AV{r) and 

En = £n + ^£n- If Eq- (9), which is the most significant piece of the present formahsm, can 
be solved analytically as in (8), then the whole problem, in Eq. (4) reduces to the following 
form 



[W„(r) + AWr,{r)f - [Wr,{r) + AV^„(r)]' = V{r) - E^, (10) 

v2m 

which is a well known treatment within the frame of supersymmetric quantum theory 
(SSQT) [14]. Thus, if the whole spectrum and corresponding eigenfunctions of the un- 
perturbed interaction potential are known, then one can easily calculate the required super- 
potential Wn{r) for any state of interest leading to direct computation of related corrections 
to the unperturbed energy and wave function. 

For the perturbation technique, we can split the given potential in Eq.(4) into two parts. 
The main part corresponds to a shape invariant potential, Eq. (8), for which the super- 
potential is known analytically and the remaining part is treated as a perturbation, Eq. 
(9). Therefore, it is obvious that SSC potential can be treated using this prescription. In 
this regard, the zeroth-ordcr term corresponds to the Coulomb potential while higher-order 
terms consitute the perturbation. However, the perturbation term in its present form cannot 
be solved exactly through Eq. (9). Thus, one should expand the functions related to the 
perturbation in terms of the perturbation parameter A, 

OO OO oo 

AV^(r;A) = ^A.V.(r), Aiy.(r; A) = ^ A,iy«(r), Ef{X) = Y.KE^, (H) 
1=1 1=1 1=1 

where i denotes the perturbation order. Substitution of the above expansions into Eq. (9) 

and equating terms with the same power of A on both sides up to 0{X^) gives 

2iy„(r)Ty«(r) - = V,{r) - E^\ (12) 

W^W(r)W^W(r) + 2Wr,{T)Wi'\T) - ^d^n\r) ^ _ ^^)^ 

V2m ar 

Wn{r)Wi'\r) + Wi'\r)Wi'\r)] - ^^^f = Vs{r) - E^\ (14) 

v2m dr 



(15) 

Hence, unlike the other perturbation theories, Eq. (9) and its expansion, Eqs. (12-15), 

give a flexibiUty for the easy calculations of the perturbative corrections to energy and 
wave functions for the nth state of interest through an appropriately chosen perturbed 
superpotential. 



III. APPLICATION TO THE SSC POTENTIAL 

Considering the recent interest in various power-law potentials in the literature, we work 
through the article within the frame of low screening parameter. In this case, the SSC or 
Yukawa potential can be expanded in power series of the screening parameter 5 as [15] 

'A\ , ^ . (A\ ~ 

1=0 

where the perturbation coefficients are given by 



^(^) = - (7) exp(-5r) = - (^) E V^.(5r)\ (16) 



VI = -1, = 1/2, = -1/6, V, = 1/24, • • • . (17) 

We now apply this approximation method to Yukawa potential with the angular momentum 
barrier 



V{r) = - ^ exp(-(5r) 



H(r) + 



+ l)/i 



+ Al^(r), (18) 



where the first piece is the shape invariant zeroth-order which is an exactly solvable piece 
corresponding to the unperturbed Coulomb potential with Vo(r) = —A/r while Ay(r) = 

A5 - {A5'^/2)r + {Ad^/^y^ - (A(5'^/24)r^ H is the perturbation term. The literature is 

rich with examples of particular solutions for such power-law potentials employed in different 
fields of physics, for recent apphcations see Refs. [16,17]. At this stage one may wonder 
why the series expansion is truncated at a lower order. This can be understood as follows. 
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It is widely appreciated that convergence is not an important or even desirable property for 
series approximations in physical problems. Specifically, a slowly convergent approximation 
which requires many terms to achieve reasonable accuracy is much less valuable than the 
divergent series which gives accurate answers in a few terms. This is clearly the case for 
the Yukawa problem [18]. However, it is worthwhile to note that the main contributions 
come from the first four terms. Thereby, the present calculations are performed upto the 
third-order involving only these additional potential terms, which suprisingly provide highly 
accurate results for small screening parameter 5. 



A. Ground State Calculations {n = 0) 

In the fight of Eq. (8), the zeroth-order calculations leading to exact solutions can be 
carried out readily by setting the ground-state superpotential and the unperturbed exact 
energy as 

(19) 

and from the literature, the corresponding normalized Coulomb bound-state wave function 
[19] 

X^n^ ir) = iViy exp [-i3r] x Ll'^' [2f3r] , (20) 



in which N^i^ 



, / / ^ is a normalized constant, /5 = , 

^n+e+i)n J (n+i+i) ^ {n+e+i)n 

and (x) = E^=o(-l)'" (n-ri)!(I,^+fc)!m! ^'" associate Laguarre polynomial function [20]. 

For the sake of calculation of corrections to the zeroth-order energy and wavefunction, 
one needs to consider the expressions leading to the first- and third-order perturbation given 
by Eqs. (12-15). Multipfication of each term in these equations by X^(r), and bearing in 
mind the superpotentials given in Eq. (7), one can obtain the straightforward expressions 
for the first-order correction to the energy and its superpotential: 



^+1 1 1 . .. , " mA 



/2m 1 



dx, (21) 



and for the second-order correction and its superpotential: 



£;(2) 
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-r^ - Wi^^ (r) l^i^) (r) 



W^^ (r 



/2m 1 



and for the third-order correction and its superpotential: 



24 



(ir, 



(22) 



Wi'^ (r 



/2m 1 



dx, 



(23) 



for any state of interest. The above expressions calculate Wi^\r), Wi^\r) and WP{r) 
explicitly from the energy corrections EI^\ E^^ and E^^^ respectively, which are in turn 
used to calculate the moderating wave function Un{r). 

Thus, through the use of Eqs. (21-23), one finds the ground state energy shift upto the 
third-order and their moderating superpotentials as 



£/0 — 



h\3N^ - L) 
4m 



5\ 



E, 



(2) _ h'N^ ibN^ -31 + 1) ^3 h'N^ ibN^ -31 + 1) 



12Am2 



16A^m^ 



5^ 



-C'O — 



h'^N^ {5N^ - 3L) {5N^ - 3L + 1) ^ h^N^ {5N^ - 3L + 1) {9N^ - 5L) 



96A2m3 

h^^N^ {5Nl - 3L + 1) {9N^ - 5L) 



-5* + 



48^3^4 



2A/2m 



8 



hNo Amr + h^NoNi 3h'^N^5-AmA 5 



(24) 



24V2m(^m)2 

where A^o = (^ + 1) ) -^i = + 2) and L = £{i+ 1). Therefore, the analytical expressions for 
the lowest energy and full radial wave function of the Yukawa potential are then given by 



En=o,e = + AS + El,'> + El'> + E^'> + ■ ■ ■ , V^„=o,^(r) ~ x'n=oAr)ur,=oAr) , (25) 



.(1) , p(2) 



.(3) 



in which 



Un=oA^) ~ exp 



2^ r fu,W 



h 



(W^^'^ (x) + WP (x)) dx 



(26) 



Hence, the exphcit form of the full wave function in (25) for the ground state is 

1 



2mA 



Am 



/+^exp(P(r)), 



(27) 



with P(r) = I]f=iPir* is a polynomial of fifth order having the following coefficients: 



Pi = 



(^+1)^(2) 



A 



Am 



9 {i + 2) 
4(£+l)' 



6 



1 

10" 



(28) 



in which d = b + j^i^^-^ig and other parameters are given in (21). 



B. Excited state calculations (n > 1) 

The calculations procedures lead to a handy recursion relations in the case of ground 
states, but becomes extremely cumbersome in the description of radial excitations when 
nodes of wavefunctions are taken into account, in particular during the higher order calcu- 
lations. Although several attempts have been made to bypass this difficulty and improve 
calculations in dealing with excited states, (cf. e.g. [21], and the references therein) within 
the frame of supersymmetric quantum mechanics (SUSYQM). 

Using Eqs. (7) and (19), the superpotential W„(r) which is related to the excited states 
can be readily calculated through Eqs. (21-23). So the first-order energy shift in the first 
excited state {n — 1) and its superpotential are 

9 



El 



(1) 



Am 



6' 



Wi'\r) ^ -^r. (29) 
2v2m 

Consequently, the use of the approximated Wi^\r) in the preceeding equation in (22) gives 
the energy correction in the second-order as 

^(2) _ h'N! (57Vf - 3L + 1) ,3 n'Nf (BNl - 3L + 1) 



We also find its supersymmetric potential 



which gives the energy shift in the third-order as 

(3) _ h'Nf (5Nf - 3L) (5Nf - 3L + 1) h'Nf (5Nf -3L + 1) (QTV^ _ 5L) 

- 9M^^ ^ + 48l%^S ^ ^^^^ 

h^^N^ {5Nl -3L + 1) (9Ar2 _ 



hNi Amr + h^NiN2 Sh'^N^S - AmA 



Therefore, the approximated energy value of the Yukawa potential corresponding to the first 
excited state is 

= +A6 + e[^^ + E^^ + ) + ■■■. (33) 

The related radial wavefunction can be expressed in an analytical form in the light of Eqs 
(21-23) and Eq.(25), if required. The appromation used in this work would not affect con- 
siderably the sensitivity of the calculations. On the other hand, it is found analytically that 
our investigations put forward an interesting hierarchy between W^^ (r) terms of different 
quantum states in the first order after circumventing the nodal difficulties elegantly, 

2v2m 

which, for instance, for the second excited state {n — 2) leads to the first-order correction 
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E. 



(1) 



h''{3Nl-L) 
4m 



r. 



(35) 



2V2m 

Thus, the use of the approximated W2^^ (r) in the preceeding equation (22) gives the energy 
shift in the second-order and its supcrpotential as 

^(2) h'Nl{5Nl-3L + l) ^, h'N^{5Nl-3L + l) ^, 

^2 = — 5 (> TTTT^ — 5 <^ ' ■ 



12^m2 



Amr + /i^iVaiVg 



6' 



24V2m(^m)2 

which leads, via Eq.(23), into the third-order energy shift 

(3) _ h^Ni{5Ni -3L){5Nl -3L + 1) , h^Nl {5Nl - 3L + 1) {9Nl - 5L) 



(36) 



E. 



h^^N^ (5A^| - 3L + 1) (9A^| - 5L) . 



(37) 



64^4^5 

where N2 = (£-1-3). Therefore, the approximated energy eigenvalue of the Yukawa potential 
corresponding to the second excited state {n — 2) is 



Er.^.,^Ei'^+AS + E^^ + Ei'^ + Ei'^ 



(38) 



Finally, from the supersymmetry, we find out the nth— state energy shifts together with 
their supersymmetric potentials as 

4m 



K = J- (> , 



2\/2m 



-r. 



^(2) ^ h\n + l + i y [5{n + l + iy-3L+l] n\n + Z + 1)^ [5(n + Z + 1)^ - 3L + 1] ^4 



nAm? 



-6-' - 



16A2m3 



5^ 



n'^{n + / + 1) \Amr + ^^(n + Z + l)(n + / + 2)1 \3h'^{n + Z + 1)2(5 - 4mAl 5^ 



2Ay/2m{Amy 
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h^jn + l + [5(n + l + 3L] [5(n + I + if - 3L + 1] . 

96A2m3 

h\n + l + [5(n + / + 1)^ - 3L + 1] [9(n + / + 1)^ - 5L] . 
h^\n + l + lf (5(n + / + 1)^ - 3L + 1) (9(n + / + 1)^ - 5L) 



(39) 



Thus, the total energy for the nth— state is 



(40) 



For the numerical results, in Tables 1-4, we list our calculated K and L-shell binding energies 
for some values of Z and compare those with the hypervirial-Pade results [8], the shifted 
large-iV expansion method [7] and the experimental values [23] for the s-statc energies Eqq 
and £^10 and also the p-state energies Eqi and En. It is observed that inspite of calculational 
simplicity, the present approach yields results as accurate as predicted by more elaborate 
hypervirial-Pade and shifted large- expansion calculation. Finally, the scope of extending 
the method to calculate oscillator strength, bound-bound transition matrix elements etc. 
which have significant importance in atomic physics is also possible. 



Table 1-3 we present our calculated iC— shall binding energies Eqq and Eqi and L— shell 
binding energies Eiq and Eu for some values of Z and compare them with the predictions of 
Lai and Madan [8] and the experimental values [17]. We quote only the Fade approximant E 
[10,11] results which provide upper bound to the energy eigenvalues. For Eiq and £^20 levels. 
We also depict our earher results [5] obtained through EWA method which provides compact 
analytic expressions only for the bound s-state energy eigenvalues. As we have used through 
out the atomic units, our energies are measured in units of 2Ry — 27.212 eV [18] is used 
.One may notice that in comparison to our earlier calculation based on EWA method, the 
present techniques gives much improved energy eigenvalues. Furthermore, our predictions 
are surprisingly close to those obtained through the use of elaborate hypervirial technique. 



IV. CONCLUDING REMARKS 
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This indicates that there is distinct advantage in using the shifted large- A?" method to similar 
calculations as it yields very accurate results yet remaining simple and straight forward. 
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TABLES 

TABLE I. Calculated K-shell energies Eqq in keV for some values of Z 



z 


EWA (Ref.5) 


Hypervirial-Pade (Ref.8) 


Shifted-A (Ref.7) 


Expt (Ref.23) 


Present Work 


3 


-0.053 34 


-0.054 15 


-0.054 14 


-0.054 75 


-0.054 056 87 


4 


-0.105 


-0.106 34 


-0.106 34 


-0.111 


-0.106 281 60 


5 


-0.178 


-0.180 08 


-0.180 07 


-0.188 


-0.180 078 08 


6 


-0.274 


-0.276 23 


-0.276 23 


-0.284 


-0.276 306 26 


7 


-0.393 


-0.395 42 


-0.395 41 


-0.402 


-0.395 579 11 


8 


-0.535 


-0.538 09 


-0.538 09 


-0.532 


-0.538 354 


9 


-0.701 


-0.704 61 


-0.704 61 


-0.685 


-0.704 983 


14 


-1.897 


-1.903 20 


-1.903 20 


-1.839 


-1.904 306 


19 


-3.716 


-3.725 45 


-3.725 45 


-3.607 


-3.727 639 


24 


-6.171 


-6.182 77 


-6.182 77 


-5.989 


-6.186 408 


29 


-9.268 


-9.282 12 


-9.282 13 


-8.979 


-9.287 593 


34 


-13.012 


-13.028 


30 


-13.028 


30 


-12.658 


-13.035 977 


39 


-17.407 


-17.424 


82 


-17.424 


82 


-17.038 


-17.435 077 


44 


-22.454 


-22.474 


38 


-22.474 


38 


-22.117 


-22.487 609 


49 


-28.157 


-28.179 


15 


-28.179 


15 


-27.940 


-28.195 740 


54 


-34.517 


-34.540 


92 


-34.540 


92 


-34.561 


-34.561 250 


59 


-41.535 


-41.566 


12 


-41.561 


17 


-41.991 


-41.585 000 


64 


-49.213 


-49.241 


18 


-49.241 


18 


-50.239 


-49.270 154 


69 


-57.553 


-57.582 


03 


-57.582 


03 


-59.390 


-57.615 917 


74 


-66.554 


-66.584 


70 


-66.584 


70 


-69.525 


-66.623 882 


79 


-76.217 


-76.250 


03 


-76.250 


03 


-80.725 


-76.294 897 


84 


-86.544 


-86.578 


78 


-86.578 


78 


-93.105 


-86.629 718 
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TABLE II. Calculated if-shell energies £^01 


in keV for some values of Z 




z 




Z 




9 


-0.012 158 


49 


-4.207 958 


14 


-0.089 499 


54 


-5.358 162 


19 


-0.282 475 


59 


-6.655 877 


24 


-0.598 417 


64 


-8.102 492 


29 


-1.044 023 


69 


-9.699 206 


34 


-1.624 349 


74 


-11.447 062 


39 


-2.343 224 


79 


-13.346 979 


44 


-3.203 631 


84 


-15.399 774 
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TABLE III. Calculated L-shell energies Eiq in keV for some values of Z 



z 


EWA (Ref.5) 


Hypervirial-Pade (Ref.8) 


Shifted- AT (Ref.7) 


Expt. (Ref.23) 


Present worl 


9 


-0.018 


-0.022 06 


-0.026 30 


-0.031 


-0.042 259 


14 


-0.116 


-0.124 92 


-0.124 92 


-0.149 


-0.130 396 


19 


-0.320 


-0.335 03 


-0.335 03 


-0.377 


-0.338 344 


24 


-0.644 


-0.665 54 


-0.665 54 


-0.695 


-0.669 125 


29 


1.096 


-1.124 48 


-1.124 48 


-1.096 


-1.128 848 


34 


-1.692 


-1.717 35 


-1.717 35 


-1.654 


-1.722 569 


39 


-2.407 


-2.448 16 


-2.450 36 


-2.373 


-2.454 212 


44 


-3.272 


-3.319 98 


-3.322 03 


-3.224 


-3.326 856 


49 


-4.281 


-4.335 27 


-4.337 19 


-4.238 


-4.342 964 


54 


-5.435 


-5.496 02 


-5.497 83 


-5.453 


-5.504 555 


59 


-6.737 


-6.803 90 


-6.805 63 


-6.835 


-6.813 316 


64 


-8.187 


-8.260 34 


-8.261 99 


-8.376 


-8.270 675 


69 


-9.787 


-9.866 54 


-9.868 12 


-10.116 


-9.877 861 


74 


-11.538 


-11.623 58 


-11.625 10 


-12.100 


-11.635 946 


79 


-13.441 


-13.532 38 


-13.533 85 


-14.353 


-13.545 871 


84 


-15.496 


-15.593 79 


-15.595 21 


-16.939 


-15.608 473 
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